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A nonequilibrium Green's functions approach to the col- 
lective response of correlated Coulomb systems at finite tem- 
perature is presented. It is shown that solving Kadanoff- 
Baym type equations of motion for the two-time correlation 
functions with the external perturbing field included allows 
to compute the plasmon spectrum with collision effects in a 
systematic and consistent way. The scheme has "built-in" 
sum rule preservation and is simpler to implement numeri- 
cally than the equivalent equilibrium approach based on the 
Bethe-Salpeter equation. 



The dynamic properties and the plasmon spectrum of 
Coulomb systems continue to attract the interest of re- 
searchers in many fields, in particular condensed matter 
theory e.g. plasma physics e.g. ||-f[| and electronic 

bilayer liquids [pi. This is due to the fact that the density 
response to an external perturbation, given e.g. by the 
dynamic structure factor S(q,u>), is a sensitive indicator 
of the state of a charged particle system which can be di- 
rectly measured in x-ray or light scattering and electron- 
energy-loss experiments. This is particularly valuable for 
strongly coupled many-body systems, such as dense plas- 
mas, metals or semiconductors at low temperature. 

In recent years there has been considerable progress 
in the theoretical account of the impact of correlations 
among the carriers on collective excitations, i.e. in 
theories which go beyond the mean-field level (time- 
dependent Hartree/ Vlasov or random phase approxima- 
tion, RPA) . It is now commonly accepted that any model 
has to obey certain consistency requirements which, in 
particular, are related to the preservation of sum rules for 
the inverse dielectric function, see Refs. [|[j4| for a discus- 
sion. Among the successful approaches, we mention at- 
tempts to construct local field corrections (see references 
in H), kinetic theory concepts to incorporate collisions 
into the dielectric function, e.g. and Green's functions 
methods originally developed by Baym and Kadanoff |uj 
and others. The latter approach is of particular interest 
as it allows for a systematic first principles treatment of 
correlations, and sum rule preservation is easily guaran- 
teed by using so-called conserving approximations for the 
Green's functions 0,|JlT). 



In most Green's functions treatments |10|,|l|-|4 11 1, lin- 
ear response theory is used which relates the density 
response function to the retarded one-particle-one-hole 
Green's function of the unperturbed system the calcula- 
tion of which is the central problem. This task is ac- 



complished by solving a Bethe-Salpeter equation (BSE), 
the quality of the results being determined by the choice 
of the four-point particle-hole-irreducible (PHI) vertex 
K. While high-level approximations for K have been 
investigated for metals at zero temperature , finite- 
temperature treatments are restricted to much simpler 
approximations Q and have to neglect vertex correc- 
tions which are important for sum rule preservation jj],[5| . 
Moreover, the collective response from a nonequilibrium 
state, which is of high interest e.g. in laser excited semi- 
conductors or laser plasmas, is completely out of reach 
in the BSE approach. 

In this letter, we present a scheme which allows i) to 
compute the linear response at finite temperature fully 
including vertex corrections, ii) the nonlinear response 
to a strong perturbation and iii) the response from an 
arbitrary nonequilibrium state. Our approach is based 
on directly computing the time-dependent density fluc- 
tuations of the electron gas under an external perturb- 
ing field, from which we obtain the density response 
function. We calculate the nonequilibrium (two-time) 
one-particle Green's functions by solving a generalized 
Kadanoff-Baym type equation with the external perturb- 
ing field included. The density fluctuation is given by the 
equal-time limit of the linearized deviation of the one- 
particle Green's function from its unperturbed value. If 
the unperturbed system is in equilibrium, sum rules are 
again preserved by adopting conserving approximations 
of the self energy jl0| . 

There is a one-to-one equivalence between levels of ap- 
proximations in the two approaches. For each choice 
of the approximate self energy in the present approach, 
its formal functional derivative with respect to the one- 
particle Green's function gives the equivalent PHI vertex 
in the Bethe-Salpeter approach. However, in investiga- 
tions involving more sophisticated approximations, our 
approach has the advantage that the self-energies are for- 
mally much simpler, and hence easier to implement nu- 
merically, than their equivalent PHI vertices. Moreover, 
our approach can be extended to compute the nonlinear 
response of a correlated electron gas. We underline that 
this efficient and conserving calculational scheme is not 
limited to the problem of plasma oscillations but is of in- 
terest also for current studies of finite temperature spin 
modes in Fermi liquids |l2| and finite temperature collec- 
tive excitations in nucleij[13[. In the following, to avoid 
confusion in terminology, we append the suffix "-BS" to 
the labels of the approximations in the Bethe-Salpeter 
approach. 
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We consider a correlated electron gas in a neutralizing 
background under the influence of an external potential 
U described by the hamiltonian H = H sys + H ex t , with 

H sys = ^ekflkflk + \ Yl T/ ('?) a k 1 +q a k 2 -q ak 2 a ki I 

k 



k!k 2 ,q#0 

H ext = ^ U(-q, t) ^ 4 +rl a k . 



q 



k 



(1) 



Here k, q are momenta, e& is the one-particle energy, and 
V(q) and U{q) are the spatial Fourier components of 
the Coulomb potential and the external potential, re- 
spectively, (a) are Heisenberg creation (annihilation) 
operators evolving with the total hamiltonian H. (Spin 
degrees of freedom are not of interest for our analysis, 
so spin indices will be suppressed.) The noncquilibrium 
properties of the inhomogeneous electron gas are defined 
by the two-time correlation functions 

G<(k + q, t i; k, t 2 ) = *<4(*a)ok+q(tx)); 

G > (k + q,ti;k,t a ) = -i(a k+q (ti)4(t 2 )), (2) 

where the statistical averaging is over the density opera- 
tor of the initial state. In particular, the density is given 
by n(q, t) — — i J2k G < (k + q, t; k, t). The time evolu- 
tion of G > and G < is governed by the Kadanoff-Baym 
equations (KBE) |h}, 

^^-e kl )G^(kii i; k 2 t 2 ) 

= ^C/(-q,i 1 )G^(k 1 -q,t 1 ;k 2 t 2 )+ (3) 
q 

^ S HF (k 1 t 1 ;kt 1 )G^(kt 1 ;k 2 t 2 ) + I*(kiti; k 2 t 2 ), 

k 

(to be supplemented with the adjoint equation), where 
E HF is the Hartree-Fock selfenergy, and the collision in- 
tegrals I< contain the short-range correlation effects (see 
below) . 

As we are interested in the dynamical response of 
the electron gas to a longitudinal electrostatic pertur- 
bation, we choose U(q,t) = L/o(i)<5q,qo- Before the on- 
set of the field, t < to, the system is homogeneous, 
G<(kxti; k 2 t 2 ) ~ <5ki,k 2 i however, for t > t , the field 
gives rise to harmonic modulations 



G^ 



G^(k 



//iq ,ti;k + // 2 q ,t 2 ), (4) 



where /ii and /i 2 are integers running from — oo to oo. 
qo enters G Mlj/i2 as a parameter and will be omitted. 
Due to the symmetry properties G< +n „ 2+n (ktit 2 ) = 



"in 



,(k+raqo,ti* 2 ) and G^ 1 , M2 (kf 1 i 2 ) 



only the matrix elements GSj, n = 0, ±1, ±2, . . . are in- 
dependent. As a result, Eq. (H) transforms into a system 

of equations for the functions G< , i.e. one of the two mo- 
mentum arguments of G^ in Eq. (||) has been replaced 



by the discrete "level" index n. Obviously, this repre- 
sentation closely resembles the multilevel/multiband ki- 
netic equations (Bloch equations) familiar from atomic 
or semiconductor optics, if written in terms of two-time 
correlation functions, e.g. [M. Only here, Goo corre- 
sponds to the spatially homogeneous state, while G< 
describes transitions of an electron from momentum state 
k + nqo at t = t\ into state k at t = i 2 , cf. definition 
(|J) . In particular, the equal-time components of Gq and 
G^q yield respectively the homogenous density compo- 
nent no(t) — —iJ2kGQ (ktt) and the field-induced fluc- 
tuations 



<5n(q, t) = Sn r 



k 



^G< (ktt). (5) 



In situations where a perturbation treatment of the ex- 
ternal field is applicable, the leading order of the Fourier 
components of the density is 8n m ~ 0(U™). Since the 
main subject of our paper is the effect of correlations 
(collisions) on the plasmon spectrum, we will focus on 
the weak-field (linear response) limit below. Then, wc 

neglect all components of G^ except G^ and G^ . Up 

to first order in the field, the equations for Gf read, for 
any fixed qo (summation over m = 0, 1 and integration 
over t from — oo to oo is implied), 



ih-^- - e k+qo ) Gf (kht 2 ) = U (h) G|(ktit 2 ) + 

E™(k* 1 )G| (ki 1 * 2 ) + 

Ef^QGloCktta) + vf m (k tl t)Gi (kit 2 ), 

~ iH W 2 ~ £k ) G ^ (kilt2) = G n( k *i*2) + 

GL(kiii 2 )£* F (ki 2 ) + 
Gf m (ki 1 t)E| (ktt 2 ) + GL(kW")s£ (kft 2 ), (6) 



whereas G^ obey the "conventional" spatially homogo- 

neous field-free equations. In Eq. (0), G^(k£it 2 ) = 

Gq (]£ + q , iii 2 ), the retarded and advanced functions 
are defined by (F denotes G or E) 

F«/£(kt 1 t 2 ) = ±e[±(t 1 -t 2 )} x 

{F> X2 (kht 2 ) - F^khh)} , Ai, 2 = 0,l > (7) 
and the Hartree-Fock selfenergy is 

E« F (ki) = S Xi V(q )J2H)Gw(ptt) 



{-i)G<(k-p,U)V(p), A = 0,1. (8) 



The selfenergies E^ and S^' 4 are of first order in 
the field and are obtained from the respective "00"- 
components by replacing one Goo a * a time by Gio and 
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summing over all terms generated this way, and En fol- 
lows from Eoo by replacing Goo by Gn. The equations 
for Goo and Gio are to be supplemented by the proper 
initial conditions for Gq corresponding to a correlated 
spatially homogeneous electron gas, see below. Further- 
more, Gio(kio^o) = 0. 

Let us now consider how the dielectric linear response 
functions can be determined from the solution of Eqs. @ 
and what their properties are. Using Eq. ([5]), we find in 
linear response, 

/>oo 

^G<(ktt) = i5n^{t) = dt x R (<lo,t,t)Uo(t), (9) 



where x R is a retarded susceptibility which, in general, 
depends on two times. If the unperturbed system is in 
a stationary state, X R {t)~£) ~ * X R (t ~ allowing to ap- 
ply the convolution theorem to Eq. (^) with the result 
(lu— dependence denotes the Fourier component), 



X 



*(qo,^) 



(10) 



which immediately yields the retarded inverse dielectric 
function and the dynamic structure factor 



(w,q ) = 1 + 



V(qo) 
1 



^G<(k, w ), 



k 



irn Q U (u}) 



ImG<(k, W ). 



(11) 



(12) 



Now, the quality of the plasmon spectrum (jl^) com- 
puted from Gio(iit2) is fully determined by the approx- 
imation for the field-free selfenergies Eoo hi Eqs. (|^). In 
particular, if Eq f = EJq = (noninteracting electrons 
gas), Eq. ( flO| ) reduces to the familiar Lindhard polariza- 
tion, \ R = n™. Further, if only the Hartree mean field 
is included, Ego = Eq , one recovers the RPA-BS result 
(full ring sum), \ R = t=VtP™ > or > equivalently, 



x R = x* + x*v x R , 



(13) 



with x* = H . Finally, with the Fock and correlation 
terms, E F and E^q, included also, one again recovers 
Eq. (|l3|), but with a more general expression for the 
proper (irreducible) polarization: 



with 



m - 



O + CO 



K 


+ 


K 




T 



(14) 



(15) 



where Eq. ([lj) starts with H R0 (first diagram) but now 
contains exchange and correlation corrections (second di- 
agram) in terms of the particle-hole T-matrix T which 
obeys the Lippmann-Schwinger equation |l4[) with the 



general interaction kernel K (see below). One read- 
ily recognizes in Eq. ([p]) the familiar field-free Bethe- 
Salpeter equation which thus is a direct consequence of 
the Kadanoff-Baym equations with weak external field 
(^). We underline that this result applies to equilibrium 
and arbitrary nonequilibrium situations (notice that all 
derivations are performed on the Keldysh contour, and 
directed lines denote full Green's functions with selfen- 
ergy insertions) |14|. 

As we demonstrate below, this connection between the 
BS and KB approaches is particularly fruitful for the di- 
electric response of a correlated electron gas: (i) as the 
BS approach is a standard formalism for the investigation 
of correlation effects, e.g. @, it can be used to classify ap- 
proximations and their properties; (ii) there exists a one 
to one correspondence between the selfenergies EJq and 
the PHI vertex K in Eq. ( |l4| ) ; (iii) based on the internal 
consistency of the Kadanoff-Baym formalism, the proper- 
ties of the plasmon spectrum are completely determined 
by the approximation for EJq: in particular, density con- 
servation of Eq (which is trivial to meet) guarantees sat- 
isfaction of the f-sum-rule . 

A valuable practical advantage of the present scheme is 
that simple approximations for Eq correspond to rather 
complex approximations for K which allows for efficient 
computation of the plasmon spectrum of correlated sys- 
tems by solving the KBE (^) . We demonstrate this below 
on the example of the (density conserving) second Born 
approximation 

E|(ktt') = iJ2 \V st (p)M (k - p,U')Gg(p,t't) (16) 
p 

/0> 



where V s t is the statically screened Coulomb potential 
(dashed lines in the diagram) and IIq-q the nonequilib- 
rium generalization of the Lindhard polarization bubble, 
n|(ktt') = -iJ2 p G%(k + p,tt')Gf Q (pt't) Proceeding 
as in Ref. the simple correlation selfenergy (|l(|), 

together with the Fock mean field, transforms into the 
following PHI vertex K in the BSE: 



K 



+ O 



□ 



(17) 



and contains contributions from particle-hole (un- 
screened) Coulomb scattering (first diagram, zig-zag line 
denotes bare Coulomb potential V), excitation of a 
particle-hole pair (second diagram) and scattering be- 
tween two particle- hole pairs (last two diagrams). Com- 
parison of Eqs. (|l6|) and ([Tt| ) reveals the familiar relation 
between E and K: K = V + <5E 00 /<5G [juM 

While it is very difficult to solve the BSE with ker- 
nel (|l7]) without further simplifying approximations, 
solving Eq. (||) with the selfenergy Eq. ( |l6| ) and the Eio 
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derived from it is quite straightforward. We performed 
numerical solutions for a strongly correlated electron gas 
in equilibrium using the numerical procedure which was 
developed before for the two-time semiconductor Bloch 
equations |16|. To create a correlated initial state, we 
run the field-free program for a time longer than the cor- 
relation time starting from an uncorrelated distribution. 
After this, the field U was turned on, where we chose a 
pulse shape for Uo(t) broad enough to cover the plasmon 
spectrum. The thus excited density fluctuation is shown 
in Fig. 1 for a 3D electron gas with Brueckner parameter 
r s = 4 and temperature fc^T = 0.69Ep (Fermi energy) 
for two wavenumbers. For comparison, also the uncor- 
related result is shown (Hartree-Fock selfenergies only 
which is equivalent to RPA-BS plus exchange). While 
5n(t) depends on the explicit form of Uo(t), obviously 
the linear response quantities e R (ui) and S(cu), Eqs. (|ll| , 
|l2| ), are independent of Uq. Fig. 2 shows the dynamic 
structure factors, corresponding to the results in Fig. 1. 
Clearly, one sees that the short-range correlations lead 
to a damping of 5n(t) in excess of the collisionless Lan- 
dau damping, cf. Fig. 1, which corresponds to a red 
shift and an additional broadening of the plasmon peak 
in the structure factor, Fig. 2. Remarkably, our numeri- 
cal scheme preserves the f-sum rule for the small (large) 
wave number to 0.03% (0.6%). In contrast, neglecting 
terms in K but keeping the selfenergy insertions in G 
leads to violation of the sum rule. For example, for in- 
clusion of the first two diagrams only (curve "1+2" in the 
inset of Fig. 2) and the first diagram only ("1") the cor- 
responding numbers are, respectively, 2.1% (0.8%) and 
1346% (416%), and for still smaller go the error increases 
rapidly ||. Clearly, the first two diagrams are the most 
important ones in producing the overall features of S{uS), 
but, as the inset of Fig. 2 shows, the last two diagrams 
may still be needed to give the correct high-energy be- 
havior. 

In summary, we have presented a new selfconsistent ap- 
proach to the dielectric properties of a correlated electron 
gas. Using an "interband" generalization and solving the 
problem in the time domain allows to take maximum 
advantage of the selfconsistent Kadanoff-Baym scheme: 
simple approximations for the collision integrals trans- 
form into complex correlation corrections in the plasmon 
spectrum with "built in" sum rule preservation. This 
scheme is straightforwardly extendable to higher order 
correlations. Moreover, it applies to arbitrary nonequilib- 
rium situations, and is easily generalized to the nonlinear 
dielectric response in strong fields. 
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FIG. 1. Density fluctuation of a strongly correlated elec- 
tron gas with for two wavenumbers. For comparison, the un- 
correlated response for one wavenumber (dotted line) and the 
exciting field (dashes) are shown too. Uf denotes the Fermi 
momentum, Ry= 13.6eF. 

FIG. 2. Dynamic structure factor (|l^) for the correlated 
electron gas of Fig. 1 (same line styles). Inset shows S for 
go = 0.62&F and contains two other approximations to the 
correlations corresponding to retaining the first diagram in 
Eq. (17) and first plus second diagrams, respectively. 
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Kwong et al. Fig. 1 
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